for i=1:6
    subplot(2,3,i)
    hold on
    scatter(bpm,c(:,i),'*k');
    plot([200 600],[0.96*mean(c(:,i)) 0.96*mean(c(:,i))],'--r')
    plot(bpm,c(:,i));
    box on
    axis([0.94*min(bpm) 1.05*max(bpm) 0.96*min(c(:,i)) 1.04*max(c(:,i))])
   
    plot([200 600],[1.04*mean(c(:,i)) 1.04*mean(c(:,i))],'--r')
   xlabel('Heart rate (bpm)');
   ylabel('Estimated T1 (ms)');
   if i==1
            legend('Measurements', 'Mean +- 4%','Location','NorthOutside')
   end
end
    